THE MAXIMUM OF BROWNIAN MOTION WITH 
PARABOLIC DRIFT 



SVANTE JANSON, GUY LOUCHARD, AND ANDERS MARTIN-LOF 

Abstract. We study the maximum of a Brownian motion with a par- 
abolic drift; this is a random variable that often occurs as a limit of the 
maximum of discrete processes whose expectations have a maximum at 
an interior point. We give series expansions and integral formulas for the 
distribution and the first two moments, together with numerical values 
to high precision. 



1. Introduction 

Let W{t) be a two-sided Brownian motion with 1^(0) = 0; i.e., {W{t))t>o 
and {W{—t))t>o are two independent standard Brownian motions. We are 
interested in the process 

W^it):=Wit)-jt'' (LI) 
for a given 7 > 0, and in particular in its maximum 

M^:= max W^{t) = max (W{t)-'yt'^). (L2) 

— oo<t<oo — oo<t<oo 

We also consider the corresponding one-sided maximum 

:= max W^{t) = max {W{t) - jt^). (L3) 

0<t<oo 0<t<oo 

Since the restrictions of W to the positive and negative half-axes are inde- 
pendent, we have the relation 

= max(iV^,iV^) (L4) 

where Ni^ is an independent copy of N.y. 

Note that (a.s.) — )• —00 as t — )• ±00, so the maxima in (1.2) and 



(1.3) exist and are finite; moreover, they are attained at unique points and 
M^,N^ > 0. It is easily seen (e.g., by Cameron-Martin) that and 
have absolutely continuous distributions. 
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1.1. Scaling. For any a > 0, W{at) = a^/'^W{t) (as processes on (—00,00)), 
and thus 

= max (W{at) - -f{atf) 

— 00<t<CXD 

= max (a^/^W{t) - a^jt^) = a^''^M./2.,, (1.5) 

— oo<t<oo ' 

and similarly 

= (1.6) 

The parameter 7 is thus just a scale parameter, and it suffices to consider a 
single choice of 7. We choose the normalization 7 = 1/2 as the standard case, 
and write M := M1/2, N := N1/2. In general, ([L5|)-((L61) with a = (27)-2/3 
yield 

= {2j)-^/^M, = {2jy^/'^N. (1.7) 

Remark 1.1. More generally, if a,b > 0, then 

max (aW(t) - ht^) = aMf,/a = (^) M. (1.8) 

-oo<t<oo^ ' ' V2o/ 

The purpose of this paper is to provide formulas for the distribution func- 
tion of M and, in particular, its moments. The main results are given in 
Section [2| with proofs and further details in Sections [3[|5| The numerical 
computations are discussed in Section [6} We use many more or less well- 
known results for Airy functions; for convenience we have collected them 
in Appendices [X1 - [B| Finally, Appendix [C] discusses an interesting integral 
equation, while Appendix [D] contains an alternative proof of an important 
formula used in our proofs. 

1.2. Background. The random variable M is studied by Barbour [3], Daniels 
and Skyrme [8j and Groeneboom [9J . It arises as a natural limit distribution 
in many different problems, and in many related problems its expectation 
EM enters in a second order term for the asymptotics of means or in im- 
proved normal approximations. For various examples and general results, 
see for example Daniels [6l[7j, Daniels and Skyrme Barbour [HIS], Smith 
[T7] . Louchard, Kenyon and Schott [14j, Steinsaltz [18], Janson [12]. As 
discussed in several of these papers, the appearance of M in these limit 
results can be explained as follows, ignoring technical conditions: Consider 
the maximum over time t of a random process Xn{t), defined on a com- 
pact interval /, for example [0, 1], such that as n — )• 00, the mean EX„(t), 
after scaling, converges to deterministic function f{t), and that the fluctu- 
ations Xn{t) — KXn{t) are of smaller order and, after a different scaling, 
converge to a gaussian process G{t). If we assume that / is continuous 
on / and has a unique maximum at a point to £ 1^ then the maximum of 
the process Xn{t) is attained close to to- Assuming that to is an interior 
point of / and that / is twice differentiable at to with /"(to) 7^ 0, we can 
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locally at to approximate / by a parabola and G{t) — G{to) by a two-sided 
Brownian motion (with some scaling), and thus max£^^(t) — ^72(^0) is ap- 
proximated by a scaling constant times the variable M, see Barbour [1]. In 
the typical case where the mean of Xn{t) is of order n and the Gaussian 
fluctuations are of order n^/^, it is easily seen that the correct scaling is 
that n^^/^(maxf X„(t) — X„(fo)) cM, for some c > 0, which for the 
mean gives Emaxt Xn{t) = nf{to) + n^/^cEM + o{n^/^), see [3 El |7|. As 
examples of applications in algorithmic and data structures analysis, this 
type of asymptotics appears in the analysis of linear lists, priority queues 
and dictionaries \13\ [Hj and in a sorting algorithm |12j . 

2. Main results 

The mean of M can be expressed as integrals involving the Airy functions 
Ai and Bi, for example as follows. 

Theorem 2.1 (Daniels and Skyrme [S]). 

2-1/3 rioc ^ 2-1/3 ydy 

EM = / = / — — (2 1) 

27Ti J_i^ Ai(z)2 27Ti Ai(%)2 ^ • ^ 

_ 2/3 Ai(t)2 + ^/3Ai(t)Bi(t) 

" io Ai(t)2 + Bi(t)2 ^'-'^ 
_ 2^/3 VsBijtf - V3Ai(t)2 ^ 2Ai(t)Bi(t) 

Jo (Ai(t)2+Bi(t)2)^ 



The expressions (2.1) and ( |2.2| (unfortunately with typos in the latter) 
are given by Daniels and Skyrme [8]. Since detailed proofs of the formulas 
are not given there, we for completeness give a complete proof in Secti on [5} 



(The proof includes a direct analytical verification of the equivalence of ( 2.1 ) 



and (2.2), which was left open in [8j.) 

By (A.l) and ( A.8| ), |Ai(iy)| increases superexponentially as y — )• ±00, 



while Ai(t) decreases superexponentially and Bi(t) increases superexponen- 



tially as t — )• 00; hence, the integrands in the integrals in Theorem 2.1 all 
decrease superexponentially and the integrals converge rapidly, so they are 
suited for numerical calculations. We obtain by numerical integration (using 
Maple), improving the numerical values in I3K5l[8l|7|, 

EM = 0.99619 30199 28363 11660 37766 . . . (2.5) 

We do not know any similar integral formulas for the second moment of 
M (or higher moments). Instead we give expressions using infinite series, 
summing over the zeros Ofc, k > 1, of the Airy function, see Appendix [A} 
Recall that these zeros all are real and negative, so we have > oi > 02 > 
see [H (10.4.94)] and Appendix |a| note that |ofc| x fc^/^, see ( |A.30[ ). 



(We use Xn ^ Vn, for two sequences of positive numbers Xn and y„, to denote 
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that < liminf„_s>oo 



Xn/Vn < limsup„_^oo Xn/yn < oo; this is also denoted 



We first introduce more notation. Let F]\f(x) be the distribution function 
of A^, i.e., Fi\f{x) := F(N < x), and let Fm{x) be the distribution function 
of M; further, let G{x) = 1-Fn{x)= F{N > x) and Gm(x) = 1-Fm(x) = 



(2.6) 



> x) be the corresponding tail probabihties. Then, by (1.4), 

Fm{x) := P(M <x)= F{N < xf = Fj^ixf. 

and, equivalently, 

Gm{x) := 1 - (1 - G(x))2 = 2G{x) - G{xf. (2.7) 

If we know G{x), we thus know the distribution of both N and M, and 
we can compute moments by 

xP'^G{x)dx, 



EMP 



P 



P 



x^ ^Gm{x) dx = p 



{2G{x) - G{xf) dx. 



(2.8) 
(2.9) 



Two formulas for the distribution function are given in the following the- 
orem. Others are given in (3.3) and Lemma 3.5 The proof is given in 
Section [3l 

Theorem 2.2. The distribution functions of M and N are Fm{x) = (1 — 
G{x)Y and Fn{x) = 1 — G{x), where 

Hi(afc) 



Gi^) = ^ E T^Ai(afc + 2^/3^), ^ > 0. 



(2.10) 



,-Ai'(a.)- 

The sum converges conditionally hut not absolutely for every x > 0. Alter- 
natively, with an absolutely convergent sum, for x > 0, 

Ai(2i/3a;) ^ 7rHi(afc) + a^^ 



G{x) 



oo 

fc=l 



-Ai(afc + 2^/\). 



(2.11) 



Ai(0) ' ^ Ai'(afc) 
The function G{x) is plotted in Figure [T] 

Remark 2.3. By (jXI]) and ( |A.30D , for any fixed x > 0, |Ai(afc + 2^/3^)| 
is usually of the order |afc|~^/^ x k~^/^, and using also ( A.31| ) and (A. 33), 
the summands in (2.10) are (typically) of the order /j-i/s-i/s-s/s _ ^-T^q 
this sum is not absolutely convergent. (For some values of k, the term may 
be smaller than because au + 2^/3^; may be close to another zero, but 
such cases are infrequent and do not prevent the series from being absolutely 
divergent.) 

'^) on the negative 



On the other hand, by (A.22), 7rHi(z) + z' 



real axis, and it follows that the terms in (2.11) are 0(/c~ ), so the series 
is absolutely convergent. Moreover, since Ai is bounded on the real axis 
(see (A.l) and (A. 4)), the sum in (2.11) converges uniformly for x > 0, 
and is thus a continuous function of x; this is no surprise since we already 
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Figure 1. G{x) 



have remarked that N has an absolutely continuous distribution, so G is 
continuous. Note also that (2.11) for x = is the trivial G(0) = 1, since 
each Ai(afc) = 0, while (2.10) does not hold for x = 0. 



The sum in (2.11 ) can be differentiated termwise and we have the following 



result, proved in Section [Sj 



Theorem 2.4. N and M have absolutely continuous distributions with in- 
finitely differentiable density functions, for x > 0, 

fN{x) - -2 ^.(0) - 2 / ^ Ai'(a,) + ^ (^-^^^ 

fM{x) = 2{l-G{x))fN{x). (2.13) 



Integral formulas for fN{x) will be given in (3.10) and (|5.10 ). The density 
functions fnix) and fuix) are plotted in Figures [2| and [sj 

Remark 2.5. In contrast, the sum 



(2.14) 



obtained by termwise differentiation of (2.10) is not convergent for any x > 
0, as will be seen in Section |3| 



Moments of M and N now can be obtained from (2.8) and ( |2.9[ ) by 
integrating (2.11) termwise. This yields the following result; see Section [4] 
for proofs as well as related integral formulas. For higher moments, see 
Remark l44l 
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Figure 2. /Ar(x) 




Figure 3. fiviix) 



We define for convenience 

ip{k) := 7rHi(afe) + 
By ([X22|), ^{k) = 0{\ak\-^) = 0{k-^/^). 



(2.15) 
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Theorem 2.6. The means and second moments of M and N are given by 
the absolutely convergent sums 

^ oo 

= ^nJJ^, - E^(^)Gi(a.) (2.16) 



fc=i 

oo 



2i/33Ai(0) 2V3 
1 vr 

3 

fc=l 



= 2i/33Ai(0) " ^ - ^(fc)Hi(a.)], (2.17) 

k=l 

E M = - , , ^ ; - ^ y r27rv?(A;)Bi(afc) - (^(A:)^] , (2.18) 

3Ai(0) 2V3Ai(0)2 2V3 Z-^ L 7 v /c; j j , v ; 

E iV' = + 2I/3 ^ v'(A;) [7ra,Gi(afc) - l] (2.19) 



k=l 



.__A^ + 2^/' E Kv'(^)Bi(a,) - a^kf] ; (2-20) 



fc=i 



2 2i/35Ai'0 4/,^ 
3Ai 0) ^ 



2 2 2Ai'(0) " 

(0) ^^vv ™fcBi(afc)--afc(/p(/c) + ^ + ^T7^ 



V9(A:) 

A:=l 

00 fe— 1 



3 ^ ai Ai(0)a^_ 

1^ (afc-aj)^ 



Numerically we have 

E{N) = 0.6955289995..., 
E(M) = 0.9961930199..., 
E(iV2) = 1.1027982645..., 
E(m2) = 1.8032957042 . . . , 

and thus 

Var(7V) = E{N^) - E{Nf = 0.6190376754 . 
Var(M) = E(m2) - E{Mf = 0.8108951713 . 



The numerical value for EM agrees with the one in (2.5). 
Further formulas for moments are given in Theorems 4.5 and 4.6 

3. Distributions 
Salminen |16[ Example 3.2] studied the hitting time 

T:=mi{t>0:x + W{t) = -I3t'^}, (3.1) 
and gave the formula |16| (3.10)], for x, /3 > 0, (with a = —/3 in his notation), 
Mt) = 2^1^ 0^1^ f; exp(2V3^2/3„^^ _ p2^3) Ai(afe + 2V^/3V3^) ^^^^^ 
fc=i ^^^> 
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for the density function of r. Note that r is a defect random variable, and 
that T = oo if and only if mint>o(x + W{t) + /3t^) > 0. By symmetry, 



oo if and only if mint>o(x + W{t) + /3t^ 
-W, and thus 



oo 



(max(H^(t) 
^ t>o 



X < 0)) = F{N,3 < x). 



Hence, choosing (3 = 1/2, 

G{x) = 1 - Fn{x) = ¥{N >x)= P(r < oo) 



fr{t)dt 



1/3 



k=l 



Ai'(afc) 



akt - ^t^) dt 



E 

fc=i 



Ai(afc + 2^/3 



Ai'(afc) 



exp(afct — ^t^) dt. 



(3.3) 



If we formally integrate termwise we obtain (2.10), from (A. 16). However, 



as seen in Remark 2.3 the sum is not absolutely convergent, so we cannot 
use e.g. Fubini's theorem, and we have to justify the termwise integration 
by a more complicated argument. 



Re mark 3.1. We have \ak\ = e{k'^/^) by ( jXsol ), and |Ai'(afc)| = e{k^/^) 
by ( fAlsIj ); further, for fixed x > 0, Ai(afc + 2^7^= 0{\ak\~^/^) = 0{k'^/^) 
by (A. 4). Hence the sums in (3.2) and (3.3) converge rapidly for each fixed 



t, because of the negative term akt in the exponent. But the convergence 
rate is small for small i, and when integrating we have the problem just 
described. 



We begin by converting the sum in (3.2) into a residue integral. Fix 
9q S (0, vr/2) and xq G (ai, 0), and let F = T{9q,xq) be the contour consisting 
of the ray {re*('^"'"^°)} for r from oo to tq := \xq\/ cosOq, the line segment 
{xQ + iy} for y G [— ro sin^Oi '^0 sin ^q] and the ray {re**^'^"^"^} for r S (ro,oo). 

For (large) integers iV G N, let Rn := (|7rA^)^/^ and let Fat := FAr(6'o,xo) 
be the closed contour obtained from F by cutting the infinite rays at r = i?7v 
and connecting them by the arc F'^y := {Rjqe^^ : 9 & [ir — 9o,Tr + ^oil- 
Note that by (A. 30), |aiv| < Rn < \aN+i\ (at least for large A^; in fact 
for all > 1). Thus, Fjv goes around the A^ first zeros of Ai; moreover, F^v 
does not come too close to any of the zeros; this is made more precise by 
the estimates in Lemma IA.2I and Lemma IA.3[ 



Lemma 3.2. Let r = be the hitting time (3.1) for /3 = 1/2 and some 
X > 0. Then the defect random variable t has the density function, for 
t > 0, 

27ri Jy Ai{z) 



frit) 



dz. 



(3.4) 
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Proof. For x > and z G F or z G F^v, Lemmas |A.1| and A. 2 yield Ai{z + 
Ai{z) = 0(1), and thus the integrand in ( |3.4[ ), ^{z) say, is bounded 

by 

$(z) = 0{exp{-2-^/^t\ Rez\ - = 0(exp(-2-^/3t| Rez\)). 

This shows both that Jp ^(z) dz is absolutely convergent, and that 



^{z)dz- <^>{z)dz-^0 asiV^oo. (3.5) 

$(z) has simple poles at the zeros of Ai, and evaluating /p^ ^{z)dz 
by residues, we see that it equals the partial sum of the N first terms of 
(3.2). Consequently, (3.2) yields J-^^ ^{z) dz — )■ frit) as N ^ oo, and the 
result follows by (3.5). □ 



axis, using the estimate in Lemma 

/3/R „ 1 



Remark 3.3. The contour F in (3^) can be deformed to the imaginary 



A.l Hence, setting z 



frit) 



(3.6) 



27r " Ai(2i/3sj) 
Moreover, this holds also for t < 0, with /r(i) = 0, since the right-hand side 



of (3.6) then easily is shown to vanish: writing it again as a line integral 
along the imaginary axis we can move the line of integration to Re 2; = a, 
for any a > ai, and for t > we may let a — )• +00, again using Lemma [A.l 

This exhibits e* ^^frit) as the i nverse Fourier transform of s 1— t- Ai{2^/^ {si+ 
x))/Ai(2^/^si). By Lemma 



A.l 



this function is integrable and in L^, and 
using ( A.l[ ) and (A.2), it is seen that so is its derivative, which implies that 
the Fourier transform e* frit) is integrable. The Fourier inversion formula 
yields 



e*'/Vr(t)e-^'*dt 



A\{2^/^{si + x)) 



Ai(2i/3si) ' 
and, more generally, by analytic continuation, 

Ai(2l/3(^ + 2;)) 



-zt 



dt 



-00 < s < 00, 



Rez> 



(3.7) 



(3.8) 



Ai(2V32) ' 

(and, in fact, for Rez > ai). This formula for the Laplace transform is 



(in a more general version) given by Groeneboom [9l Theorem 2.1], where 
e* ^^frit) is denoted /ii/2,a;(*); see also (5.5) below. Conversely, this formula 
from [9j yields by Fourier inversion (3.6) and (3.4), so we could have used 



it instead of (3.2) from [16j as our starting point. We give an alternative 
proof of (3.8) in Appendix [D| which thus gives us a self-contained proof 

and our Appendix |D| use 



of Lemma 3.2 (Groeneboom [9j, Salminen 
similar methods. See also Appendix O for another approach.) 



Remark 3.4. For our purposes we consider only x > in (3.1). For x < 0, 



the hitting time is a.s. finite; its distribution is found in Martin-Lof [15] 
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Lemma 3.5. For x > 0, 



G(x) = - 
^ ' 2i 



1 f m{z)A[{z + 



Ai{z) 



dz. 



(3.9) 



Proof. We have G{x) = P(r < oo) = frit) dt. Now integrate with 
respect to t and interchange the order of integration, which is allowed by 
Fubini's theorem and the estimate Lemma [A. 1[ which implies that for z £T, 
the integrand in (3.4) is bounded by 0(exp(— cx|2;| 

_ t^/Q)) . The result 

follows by (A. 16) (and a change of variables t = 



□ 



The integral in (3.9) is absolutely convergent and converges rapidly for 



any fixed x > by Lemma A.l and (A. 21). We denote the integrand in 
by ^-(2) = '^{z;x). 



XT] and [X2 

7^1/^)) 



Proof of Theorem \2.^ Fix x > 0. By (A. 21) and Lemmas 
for z = i2ive*('^+'')Vith \(p\ < f, ^(z) = 0{\z\~'^ exp{-x\z\ 
follows that Jp^ ^{z) dz - Jp ^(z) dz as 00, and thus G{z 

^ limAr_^oo /p^^ ^{z) dz. Evaluating Jp by residues, we obtain (2.10). 



It 



To obtain (|2.11|), we take out the first term in the expansion (|A.20|) of 
Hi(z) 



G{x) 



and write (|3.9|) as 
1 



(Hi(z) 



+ 7r"iz"i)Ai(z ■ 



2^/'^x) 



dz 



dz. 



1 f Ai{z + 2^/^x) 
/r Ai(z) 2m zAi(z) 

The first integral can be converted to a sum of residues by the argument 



just given for (3.9), which yields the sum in (2.11). Indeed, we have better 



estimates now, and the resulting sum is absolutely convergent as seen in 
Remark 12. 3[ 

For the second integral we instead close the contour on the right, by a large 
circular arc {i?e**} for t from vr — 
that the error tends to as i? - 
no zeros, so the only pole is at z 



60 to — (vr — ^o); it follows by Lemma A.l 



00. Inside this closed contour, Ai has 
where the residue is Ai(2i/3x)/Ai(0). 
The result follows, noting that we go around this contour in the negative 
direction. □ 



Remark 3.6. We may also use the expansion (A.20) of Hi with more terms. 



In general, subtracting the sum with L terms in (A.20) from Hi in (3.9) 



yields an integral that can be converted to a sum of residues as above; this 
sum is similar to the ones in (2.10) and ( 2.11| ), and the terms are now of 
order k~^~^^. We also have the integral with the subtracted terms; this 
is a linear combination of terms of the type J-^z~^~^A[{z + 2^^^x)/ A[{z), 
which as above equals — 27ri times the residue at 0, so this integral can be 
written as a combination of derivatives of Ai at 2^/^x and 0; by the equation 
Ai"(z) = zAi(2;), and successive derivations of this equation, the result can 
be written as pi{x)Ai{2^^^x) + p2{x)A[' {2^^^x) for some polynomials pi and 
P2 (depending on L), whose coefficients are rational functions in Ai(0) and 
Ai'(O). We leave the details to the reader. 
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Proof of Theorem 2.4 If | arg z\ < vr— 5 and \ z\ > 1, say, then by Lemma A. 3 
|Ai'(z)/Ai(z)| X l^p^. More generally, using also Ai"(2:) = zAi{z), and dif- 
ferentiating this equation further, by induction, Ai^"^\z) / Ai{z) = 0(|z|"^/^) 
for every fixed m > 0. 



It follows by Lemma A.l and (A. 21) that for every fixed m > 0, x in a 



fixed interval {xq,xi) with < xq < xi, and z gT, 



dx^ 



;'if{z;x) = O 



^|m/2g-ca;|2:|i/2 



o 



^\m/2-l^-cxo\z\^^^ 



Consequently, we can differentiate (3.9) under the integral sign an arbitrary 



number of times; this shows that G is infinitely differentiable on (0, oo) and 
that has an infinitely differentiable density /at = —G' given by 

r Hi(z)Ai'(z + 2^/^x) 



fNix) 



-G'(x) 



Ai{z) 



dz. 



(3.10) 



This integral can be evaluated as a sum of residues as for (2.11) in the 
proof of Theorem 2.2 above, by first adding 7r~^z~^ to Hi(z), which yields 



(2.12). Alternatively, and perhaps simpler, (A. 5) and (A. 31) imply that. 



uniformly for < x < xq for any fixed xq > 0, Ai'(afc-|-2-'^/^x) = 0{\ak\^/'^) 
0{Ai'{ak)), and thus the terms in the sum in (2.12 ) are 0(|7rHi(afc) +a^"'^|) 

-4 



Oi\ak\-^) = 0(/t-8/3). Hence we can integrate the sum in (12.12^ termwise; 
equivalently, we can differentiate (2.11) termwise, which yields ( 2.12[ ). 
The result for M and follow from Fm(x) = F/v(x)2 = (1 - G{x)f. 

□ 



To see that the sum (2.14) does not converge for any x > 0, let y := 2^^^x. 
Take x = |afc| — y in (A.5). Since then, by Taylor's formula and (A.30), 

7r(4A; - : 



y\ak 



1/2 



+ 1 



l{\a,\-yf/' = 
we obtain 

Ai'(afc + y) 
and thus 

Ai'(afe + y) cos(7rA; - y(37r/c/2)^/3) _^ 



y(37rA:/2)i/3 + o(l), 



-vr 



i/>fe|^/^(cos(7rA; - y(37rA;/2)i/3 + o(l)) + o(l)) 



cos 



(y(3^fc/2)i/3)+o(i) 



Ai'(afc) cos(7rA:) + o(l) 

Let In be the interval [c(27rn)3, c(27rn + 1)^], with c := 2y^^/{3TT). For 
k £ In, we have y(37rA;/2)^/3 ^ [27rn, 27rn + l] and thus Ai'(afc + y)/Ai'(afc) > 
cosl + o(l) > 0.5, if n is large. Further, by ( |A.19[ ) and (|A.30|), Hi(afc) 



TT 



TT -^(37r/2) ^/3^ 2/3_ Hence the term in (2.14), tk say, satifies, 



for some constants ci, C2 > and k £ In, 

tk > ciA;"^/^ > 



C2n 



Since there are 0(n ) integers in /„, the sum over them is 6(1), and thus 



the sum in (2.14) diverges. (The case x = is simple.) 
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4. Moments 

Lemma 4.1. For every fixed p > 1, uniformly in a >0, 

/•oo 

/ xP-^\Ai{x-a)\dx = OiaP~^/^ + l), (4.1) 
Jo 

/ xP'^\Ai{x-a)\'^dx = OiaP~^/^ + 1). (4.2) 
Jo 

Proof. For < x < a we have xP~^ < a^^^ and | Ai(x — a)\ = 0{{a — x)^^/"^) 
by (ig, so 

£xP-^\Ai{x - a)\ dx = (^aP-^ j\a - x)-^'^ dx^ = 0{aP-^a^/^), 

and similarly xP~-^|Ai(x — a)p dx = O [dP-^a^l'^). For lar ger X we use the 
rapid decrease in (A.l), which implies 



x^'~^|Ai(x - a)| dx 



/•oo 

/ {y + a)P-^\A\{y)\dy = 0{l + aP-^) 
Jo 



and similarly for xP -'^|Ai(x — a)p dx. The result follows. 



Proof of Theorem\2.(\ Write (2.11) as 



G{2'^/^x) = c{k)Ai{x + Ofc), 



□ 



(4.3) 



k=0 

where we for convenience define = and 
c(0) := ^ 



Ai(0) ' 

Ai'(afc) Ai'(ofc) 



By (A^), (A.31) and (A.30), 



\c{k)\=0{\a,\-^-'/^)=0{k-^y'), k>l. 



(4.4) 



C>(|ofcp/4) = 0{k^/^), k>l, and thus 

the sum in (|4.3|) converges absolutely in L-'^((0, oo), dx), so it may be inte- 



By (4.1), ||Ai(x + afc)||i 



l((0,oo),dx) 



grated termwise. Consequently, using (2.8) and (A.23), 

POO j-OO °° 

EN= G(x)dx = 2^1/3 / G'(2-i/3x)dx = 2-1/3 Vc(A:)AI( 

^0 Jo ,„_n 



We have AI(0) = J^^ Ai(x) dx = 1/3 [1, 10.4.82], see (|X23]), a nd, fo r k>l, 
Al(afc) = -7r Ai^(a fc)Gi (afc) b y ^A^ since Ai(afe) = 0. Thus ( pie] ) follows, 
and so does by dXlTl). 
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Lemma 4.1 and (4.4) similarly imply that (4.3) converges absolutely also 
in L^((0,oo), dx). Hence, the integral G{x)^ dx can be obtained by 
termwise integration in (|4^, using ([RS]), ( |B.16[ ), ( |B.20D , ( |B.24| ): 



Gix f dx = 2-1/3 / ^(2-1/3^)2 
Jo 

= 2-1/3 ^^c(A;)c(^) y Ai(x + afc)Ai(x + a^) dx 



fc=0 1=0 



2-i/3c(0)2(Ai'(0))2 - 2-1/32 ^c(0)c(A;) 



Ai(0)Ai'(afc) 



k=l 



+ 2-1/3 ^c(A;)2(Ai'(afc))^ 



fc=i 



Ui(o)j h 



Thus, by (2.9) and (2.17), 



/"CXD /"OO /"OO 

EM = 2 G{x)dx- I G{xfdx = 2EN- G{xf dx 
Jo Jo Jo 



J!^ _ 2-1/3 fMo)y_^f^ 

3Ai(0) ^ VAi(0)y 21/3 A. 



27rBi(afc) - 27rHi(afc) h (p{k) 



and (2.18) follows by the definition of (p{k). 

Similarly, (4.3) converges absolutely in Li((0, 00), x dx) and L^((0, 00), x dx) 
too, and termwise integration in (4.3) yields, using (B.15), (A. 26), (B.6), 

(|b:25|, dBiQi), miT^, 



roo f-QO 

/ G(x)xdx = 2-2/3 / G(2-i/3x)xdx 
Jo Jo 

/"OO 

= 2-2/3 V / c(A;)Ai(x + ak)x dx 



fc=o ■ 

00 



2-2/3 ^ c(A;) (-Ai'(afc) - afcAI(afc)) 



fc=0 



- ^,^ + 2-2/3 J^^(fc)(-l + .a.Gi(a.)); 



k=l 



f'OO f'OO 

/ G(x)2xdx = 2-2/3 / G(2-i/3x)2xdx 



2-2/3 ^^c(A;)c(j) / xAi(x + afc)Ai(x + Oj) dx 
fc=o j=o ''^ 
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22/33Ai(0) 
1 



+ 



22/3 



k=l 



2Ai^(0) 
Ai(0)a| 



22/3 



k=l 



k=l j=l 



a,)2 



By ([2^, EN^ = 2 G{x)x dx, and (|2l9])-((23oJ) follow. Similarly, EM"^ = 
4 G{x)x dx-2 G{xfx dx = 2EN'^-2j^ Gix)^x dx, and d^IMj ) fol- 
lows. The numerical evaluation is done by Maple, using the method dis- 
cussed in Section |6l □ 

Remark 4.2. The formula above for G{x)'^ dx may be simplified. In 
fact, (Ai'(0)/Ai(0))2 = Efcli ^ by ( |B.27D , and thus the formula can be 
written 







dx = 2-1/3 J^iv^ik) - 1 /akf = 2-1/3^2 ^ jji(afc) 



(4.5) 



k=l 



k=l 



This can be seen as an instance of Parseval's formula, see Remark B.l 



However, although simpler than our expression above, this sum converges 
more slowly and is less suitable for our purposes. 

Remark 4.3. As is well-known, see [H 10.4.4-5], Ai(0) = 3-2/3/r(2/3) = 
3-i/6r(i/3)/(2vr) and Ai'(O) = -3-V3/r(l/3) = -3V6r(2/3)/(27r). We 
prefer to keep Ai(0) and Ai'(O) in our formulas. 

Remark 4.4. Higher moments can be computed by the same method, with 
Airy integrals evaluated as shown in Append ix [B| but in order to get con- 
vergence, one may have to use a version of (2.11) with more terms taken 
out of the expansion (A.20) of Hi, as discussed in Remark 3.6 We do not 
pursue the details. 



We can also give integral formulas based on Lemma 3.5 
Theorem 4.5. The moments of M and N are given by, for any real p > 0, 



ENP 



EMP 



-p2 



'P/3-l, 



+ p2 



-p2 



-p/3-2 



-p/3, 



xP'^Ai(z + x) dx dz 
Ai(z) 

xP'^Ai{z + x) dx dz 



r Jo 



r Jo 



00 



Ai{z) 

'^Ai{z + x)Ai(u; -|- x) dx 



r Jr Jo 



Hi(z) m{w) 
Ai(z) Ai(u;) 



dz dw. 



Proof. Immediate from (2.8)-(2.9) and Lemma |3.5| (with a change of vari- 
ables X — >• 2-i/^x). The double and triple integrals converge absolutely by 
Lemma rO and ^Al^. □ 



For integer p, the integrals over x in Theorem 4.5 can be evaluated by the 
formulas in Appendixpl In particular, by (B.l) and (B.22), 
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EM 



2_i/3^ r Ai{z)m{z) 



+ 2 



-7/3 



Ai{z 
Ai{z)Ai'{w 



dz 



Ai'{z)Ai{w) m{z)m{w) 



r Jr 



w 



Ai{z)Ai{w] 



dzdw. (4.6) 



Although there is no singularity when z = w in the double integral in (4.6), 



it may be advantageous to use different, disjoint, contours for z and w. 
Remember that T = T{9q,xo). We choose 9i E {9o,7r/2) and xi £ (xo,0), 
and let V := T{9i,xi). Then F and V are disjoint; moreover, if z E F and 
w GT' , then 

Iz — > cmax(|2;|, Iwl) (4.7) 



for some c > 0. Furthermore, (4.7) holds also if z E F„ and w E F^y_with 

and 



4.5 



M > 2N. We can replace the double integrals J-p Jp in Theorem 
(4.6) by /p/p/. Taking residues, this leads to another formula with sums 
over Airy zeros. 



Theorem 4.6. 



EN = 2-1/3^2 ^Hi(afc)(Hi(afc) - Bi(afc)), 



k=l 



EM = 2~^/V2^Hi(afc)(Hi(afe) - 2Bi(afc)). 



(4.8) 



(4.9) 



k=l 



These formulas are closely related to (2.17)-(2.18). They are simpler. 



but less suitable for numerical calculations since they do not even converge 
absolutely; the terms in the sums decrease as k'^/^ by (|A32J) and (|X33]). 
(However, they alternate in sign, and the sums converge.) The formulas 
([48]) and (|49]) are what we obtain if we substitute ( |2.10| ) in ( [2^ and ([2?9| 
(with p = 1) and integrate termwise; however, since the resulting sums are 
not absolutely convergent, termwise integration has to be justified carefully, 
and we use a detour via complex integration. 



Proof. Let 



/•oo 

Q{z,w) := / Ai{z + x)Ai{w + x) dx. 
Jo 



The integral converges absolutely by ( A.l ) for all complex z and w, uniformly 



in compact sets, and thus Q is an entire function of two variables; moreover. 



(B.5) and (B.9) yield the explicit formulas 

Q{z,z) = Ai'{zf-zAi{zf, 

Ai{z)Ai'{w) - Ai'{z)Ai{w) 
z — w 



Q{z,w) 



z ^ w. 



(4.10) 
(4.11) 



By Theorem |4.5 | and (4.6), using F' as discussed above, 
_2-4/3- / 



EiV 



Ai(z) 
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AI(z)Hi(z) 
Ai(^) 



r Jv 



Q(z,w)Hi(z)Hi(u;) 
Ai(z)Ai(w) 



First consider the simple integral, Jp ^(-z) dz say. It follows from Lemma A. 3 
and (A. 21) that Jp ^(.z) dz — Jp "^(-z) — s- as n — t- oo, and we find by the 
residue theorem applied to r„, letting n 



(A.17), 



AI(z)Hi(z) 
Ai(^) 



oo, together with (A. 26) and 
AI(afc)Hi(afc) 



k=l 



Ai'(afc) 
27r2z^Gi(afc)Hi(afe) 



k=l 

oo oo 

= -iTT^i Bi(afc)Hi(afc) + 27r^i ^ Hi(afc)2. 
fe=i fc=i 
The sums converge, e.g. by the argument just given, but only the final sum 
^Hi(afc)^ converges absolutely, by (A.32)-(A.34). This yields the result 
([48]) for N. 

Next consider the double integral. If z G F and u; G F', or z G F„ and 

[th m > 

Q{z,w) 



w G F^ with m > 2n, then by (4.7), (4.11) and Lemma A. 3 

Ai'(u;) 



M{z)Ai{w) 



< 



< 



C 

\z\ + \w\ 
C 



Ai'(z) 



Ai(z) 



+ 



\z\^'^ + \w\ 



\z\ + \w 

It follows that /p Jp, — /p Jp, 
theorem for first L and then L, , with 

Jin -"i m 

integral above equals 

n 2n 



Ai{w 



'1/4,^1-1/4^ 



as m > 2n — 7- OO. Using the residue 
m = 2n, we find that the double 

\ - \ - Qjaj, afc)Hi(aj)Hi(afc) 



j=i k=i 



M{aj)M{ak) 



By ( |4.11D , Q(aj,afc) = when j 7^ /c, and Q{ak,ak) = Ai'{akf by ( |4.10D . 
Hence the double integral equals —Air'^ YlT=i ^K^'k)'^- 

The result (4.9) follows by combining the terms. □ 

5. Proof of Theorem 12.11 
Proof of (|2.1[). We shall use Groeneboom [9J. Fix 7 > 0. (We may choose 



e.g. 7 = 1/2 as in other parts of this paper by (1.7), but for ease of compar- 
ison with [9j, and because we find it instructive to see how the homogeneity 
works, we write the proof for a general 7.) Fix also s > and define 



V, 



Thus, Vo = Ny, and / 



max{W{t) - -if 

t> — s 



(5.1) 



My as s 



00. 
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For t > —s, the process Ws{t) := W(t) — W{—s) is a Brownian motion, 
starting at at time —s. Define 

M = Ms:= max(W,(t) - 7*^) = Vs - W{-s), (5.2) 

and let r = be the (a.s. unique) time with Ws{t) — = M. Note 
that M > Ws j—s) — 75^ = —75^ and t > —s (strict inequalities hold a.s.) 



Groeneboom [HI Corollary 3.1] applies to Ws{t) — 7t^ (with s replaced by —s 
and X = —75^), and shows that r and M have a joint density, for t > —s 
and y > X = —75^, 

/r,M(*'y) = exp(-i7^(i^ + s^) + 27s(y + ^s^))h^^y+^,2{t + s)k^{t) 

= exp(|7^s^ + 2jsy)h^^y_^^s^{t + s)g^{t). (5.3) 

where the functions h^^y_^_^g2, k-y and are given in p]. Integrating over 
t > — s we find, with h^^aif) = for t < and §-^{1) := g.y{—t), the density 

fM °f ^ fo'^ y > ~'>'*^' 

f{t,y)dt = exp{'l-f^s^ + 2-fsy) / /i^ ,^+^^2(i + s)c/^(t) dt 

-s ^ —s 

= exp{^j^s^ + 27sy)/i^^ 2^4.^32 * ^^(s). (5.4) 

By IS Theorem 2.1], for a > 0, /i-va > and /i^^a has the Laplace 
transform, for A > 0, (see also Remark |3.3[ ) 

« := (27=)-'"A^ (5.5) 



00 







Letting A \ we see that 

/i^,a('u)dn = Ai((47)^/^a)/Ai(0) < 00, 



so h^^a £ L-'^(M) and (5.5) holds for all complex A with Re A > by analytic 



continuation. In particular, h^^a has the Fourier transform, see ( |3.7[ ), 
■ — , ^ ■ , , Ai((47)V3a + i(272)-i/3 ) 



Furthermore, by (5.3), > 0, and by Corollary 3.1], g^ G L^(M) with 
the Fourier transform 

foo 2V3^-i/3 

Hence, for y > —75^, h^^y_^_^g2 * has the Fourier transform 
^ , > ,1/3 _i/3 Ai((47)'/»fa + T»') + i(27')-'-'M 
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Since |Ai(i2/)| x \y\ exp(^|y|^/^) as y — t- ±oo by (A. 3), (5.6) implies 
|^(a;)| X exp(— 3^|a;p/^) ; consequently^ G L^(M). Furthermore, 

hj^ai^) is bounded for each a (by /i^^a G or by Lemma A.l ), and thus the 
product h^ y^^g2{uj)g~^{—uj) G ^-'^(M). Consequently, the Fourier inversion 
formula applies to (|5.7|), and (|5.4|) thus yields, for y > —75^, 



2I/3 -1/3 



, Ai((47) V3(y + ^s2) + i(2^2)-l/3^) 



d(^ 



(47) 



1/3 



27r 



exp(|7^s^ + 27s?/) / e 



Ai(i(272)-i/3^)2 

.(27^)V3, Ai((47)V3(^^^^2)^ -^^ 



Ai{iv) 



dv. 



(5.8) 



where the double integral is absolutely convergent by Lemma |A.l 



Ee 



Multiplying by e^^ and integrating, we obtain for any z € C the following, 
ale integral is 

e''fM{y)dy 







(47)^ 



1/3 roo foo 



2-K 



1 

2^ 
1 

2^ 



x=0 ^ v=—oo 
00 



,,(27^)i/., Ai((47)^/=^x + it;) 

Ai(it>)2 



du dx 



g-|72s3-7s22 



g-|72s3-7s22 



^(2+27s)(47)-l/3x+i.(27^)l/3^ Ai(x + iv) 



v=—oo J x=0 
00 is{2-y'^y^^v foo 



A\{iv) 



x=0 



Ai(w)2 

(2+27s)(47)-l/3^A.^^^ .^^ dxd?;. 

(5.9) 



Since Ai is bounded on the negative real axis by (A.4), Lemma B.2 implies 
that, for Re z > 0, 



/ e^*Ai(t) dt = e*'/3 - / e^*Ai(t) dt = e*'/^ + o[ 



e^"^* dt 



Moreover, (A. 3) implies that for y G M and 2: G C, 

e''Ai{w)dw = o[e'^\y\''"+\'^'\\y\ 



Hence, if Re 2; > (47) say, and f G M, then, using Cauchy's integral 
formula on a large rectangle with vertices {0, iy, R,R + iy} and letting R — )• 
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oo, using (A.l) to control the tails, 

foo+iv 



e''^Ai{x + iv) dx 







ivz 



e'"'A[{w) dw 

f'oo f'iv 

/ e''"Ai{w) dw - e-™^ / e^"'Ai( 
Jo Jo 



w) dw 



Hence, if Hez + 27s > 1, (5.9) yields, using (A.l) again for the error term, 
1 



Ee 



zM 



27r 



-e 3 



^■y'^s^-js^z 



00 



,{2+27s)3/(127)~jt)(z+27s)(47)-i/3 



00 Ai{ivf 

00 



dv 



1 ^,5^ 
0127""" 2 

27r 



-j(47)- 



00 
1/3, 



|Ai(it;) 



^^|«|3/2 + |Im2(47)-l/3|H 



+ Q/g-§72s3-7s2Re2+0(|Im2|2) 



) 



, Ai(ii;)2 

In particular, if 27s > 2, say, this holds uniformly for \z\ < 1, and we may 
differentiate the analytic functions at z = and obtain 



EM 



1 

2^ 



"i']^"f" d. + 0(e-^^-^+^-^). 
Ai(?i;)2 ^ ^ 



Since E = E(M + W^(-s)) = E M by ([5^, we find by letting s ^ 00 and 
choosing 7 = 1/2, 



E M = hm E Vs 



2-1/3 /-oo 

27ri 



Ai(it;) 



df , 



which is (2.1 ). 



□ 



Remark 5.1. Setting s = and 7 = 1/2 in (5.8), we obtain another formula 
for the density of N: For y > 0, 



fN{y) 



2I/3 rco Ai{2^/^y + iv) 



2n 



Ai{iv) 



dv. 



(5.10) 



By residue calculus, as with similar integrals in e.g. the proof of Theorem 2.2 



this may be written as a sum of residues of 2^/^Ai(2^/^y + 2;)/Ai(z)^ at the 
poles ak', however, now the poles are double and we omit the details. 



The integral formulas (2.1)-(2.4) can be transformed into each other by 
properties of the Airy functions as follows. 



Proof of (2.4). The integrand in (2.1) is analytic except at the zeros of Ai, 
which lie on the negative real axis. Furthermore, by (A.l), |Ai(2:)| is expo- 
nentially large, so the integrand in (2.1) is exponentially small, as \z\ — )• 00 
with 7r/3 + 5 < I arg(z)| < vr — (5; in particular when 7r/2 < | arg(z)| < 2tt/3. 
Hence, we can deform the integration path from the imaginary axis to any 
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reasonable path in this domain. We choose to integrate along the two rays 
from the origin with argz = ib27r/3, and obtain thus 



EM 



2-1/3 
27ri 



zdz 

-27ri/3oo Ai(z)^ 



2-1/3 i-oo 



2m 



which by the formula [1, 10.4.9] 



Ai(e±2^V3t)2 ' 

(5.11) 



Ai(ze^2W3) = le^"*/3(Ai(2) ^iBi(z)) 



(5.12) 



yields 

EM 



2-1/3 i-oo 
22/3 



/Hdt 



{ki{t)^iB\{t)Y 



22/3 
22/3 

vr 



g±2^i/3Ui(i)±^Bi(t))^ 

(Ai(t)2 + Bi(t)2)2 
^ (-1 ± ^/3i) (Ai(t)2 - Bi(t)2 ± 2iAi(t)Bi(t)) 



± 2(Ai(t)2 + Bi(t)2) 

°° ^/3Ai(t)2 - V3Bi(t)2 _ 2Ai(t)Bi(t) 



(Ai(t)2 + Bi(t)2)^ 



tdt, 



which proves (2.4). 



□ 



Proof of (2.2). To prove (2.2), we as above deform the integration path in 

m the 
lite in1 

dz Bi(z) 



(2.1) and integrate along the two rays from the origin with argz = zb27r/3 
and obtain (5.11). We now use the indefinite integral 



Ai(.)2 - "Ai(^ (^-^3) 

given by [3] (and easily verified by derivation, using the Wronskian Ai(z)Bi'(z)- 
Ai'(z)Bi(z) = l/vr [T| 10.4.10]). We have by and [1, 10.4.6] 



Bi(ze 



27ri/3\ 



e"'/^Ai(ze-2-V3) + g- 



1 -7ri/6 



(3Ai(z) + iBi(z)) 



and thus, by (5.12) again, 
Bi(ze2'^*/3) 



Ai(ze2'^V3) 



,/2 3Ai(z) + iBi(z) 
M{z) - iBi(z) 



vri/: 



4i 



Ai(z) 



Ai(2;) - iBi(z)' 



By (A.l) and (A.8), this converges rapidly to z as 2; — )• oo along the positive 
real axis. Hence an integration by parts yields 

2W3«, _^ /Bi(z) N ie2-/3oo '■^""'Z' 

/ . . , ,0 dz = vrzl 



Ai(z) 



L VAi(z 

+ e^^'l^ 



i I dz 







~ 'Bi(z) 
vrl 



Ai(z) 



i I dz 



vr 



4i 



Ai(t) 



Ai(t) - iBi(t) 



dt. 



(5.14) 
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The integral along the line from e~'^'^'^^^oo to equals —1 times the complex 



conjugate of the integral in (5.15), so we obtain from (5.11) and (5.15), 



EM 



2-1/3 



Im 
71" .10 



z dz 
Ai(z)2 



^2/3 r + ^/3Ai(t)Bi(0 







Ai{ty + Bi{ty 



dt, 



which is (2.2) 



□ 



Proof of (2.3). Follows similarly by (5.11) and (5.14), we omit the details. 

□ 

6. Numerical computation 

The sums in Theorem |2.6| converge, but rather slowly. For example, in 
([2l6l), \ip{k)Gi{k)\ = \ip{k)\ \B[{k) - Hi(fc)| ~ ck'^^^ for some c > 0, see 



the asymptotic expansions below. 



To obtain numerical values with high accuracy of the sums in (2.16) 



(2.21 ), we therefore use asymptotic expansions of the summands. More pre- 
cisely, for each sum, we first compute X]fc=i numerically (All computations 
are done by Maple.) We then evaluate YI200 replacing the summands 
by the first terms in the following asymptotic expansions. (We provide here 
only one or two terms in each asymptotic expansions, but we use at least 
5 terms in our numerical computations.) We note that Bi(afc) alternates 
in sign, so for terms containing it, we group the terms with k = 2j and 
k = 2j + 1 together, for j > 100. The resulting infinite sums are readily 
computed numerically. (They can be expressed in the Riemann zeta function 
at some points.) 

We use the expansions, see [H (10.4.94), (10.4.105)], ( |X20l ), ( [2l5| ), ^A^ , 



32/^2/3 

22/3 
1 _1 



k-'/' + ... 



-a. 



Hi(ofc) 

Bi(afc)~(-l)^- 



22/3 



TT 



-2al 



IT 



32/37^5/3 



(A:-2/3 + V5/3 + j 



+ ••• 

2I/6 



2II/3 



38/37^8/3 
:k~'/' + ..., 



+ 



31/6^2/3 

51 (i) := ip{2j)Bi{a2j) + ^{2j + l)Bi(a2j+i: 



92(i) := a2jV3(2j)Bi(o2j) + a2j+w{2j + l)Bi(a2j+i) ~ 



17 • 3^/6 
162 7rio/3 

13 • 3^/6 



r'"' + . . . , 



162 7r8/3 



-19/6 



+ . . 



Thus, for example, EA^ 



qN 
•^0 



+ 5f + 5. 



2 , where 



0.6955290109 . 
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5f 



vr V ^ NT-.- / N 



21/3 



fc=200 

oo 



21/3 



^ 51 (j) = 0.5317... •10-^ 



i=100 



2 ■ 21/3 



^ ^ (^(fc)Hi(afc) = -0.16722.... 10-^ 



A;=200 



yielding E(iV) = 0.6955289995 . . .. Similarly, EM = S^^ + S'f + 5f , with 



5, 



M 



0.9961930092 . . , 



5f 



vr 
'21/3 



^ 2v9(fc)Bi(afc) = 25^ = 0.10635 . . . • 10"^ 



fc=200 

00 



2 = 0.20... -10-^^, 



fc=200 



yielding EM = 0.9961930199 . . ., which fits with ([2^. Actually, we have 
15 digits of precision with the expansions we have used. 



For the second moments, we compute the sums in (2.20) and (2.21 ) in the 
same way. For the double sum in (2.21 ), we compute the sum with k < 199 
exactly, and find using asymptotics the tail sum 

27/3 y = 0.4.... 10-10. 

(We may, for example, use the first term asymptotics for (p{j) and aj for 
j > 200, since the sum is small and only a low relative precision is needed.) 

The (complementary) distribution function G{x) and the density function 
/Ar(x) may, for any given x > 0, be computed to high precision from ( |2.11 ) 
and (2.12) in the same way; for the tails of the sums ^^^^200 + 
Ai' (ak) and Yl'kL20o'Pi^)-^^'i'^k + '^^^^x)/Ai'{ak) we use the asymp- 
totic expansion oiip{k) given above together with, see [U (10.4.94), (10.4.97), 
(10.4.105)] and (|X4])-(|X5]), 

Ai(afc + 2^/3^) ~ (-l)'=+i-^^J^sin((3vr)i/3^A;i/3) fc-^/^ ^ _ ^ ^ 



Ai'(afc + 2^/3^) ~ (-1 



\k+l 



31/6^2/3 
3I/6 



2i/6vri/3 

for fixed x, and, as a special case, 



cos 



{{37r)'/''xk'/'') k'/^ + ... 



Ai'(afe)~(-l)'=+i- 



3I/6 



2i/6vri/3 



/fcVe + . . . 



The distribution and density functions of M then are given by ( |2.6| )-(2.7) 
and ( [27L31 ). 

As an illustration, we plot the tail sum (from k = 200) for /Ar(x) in 
Figure |4j 
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Figure 4. Tail sum from k = 200 for /Ar(x) 

Appendix A. Some Airy function estimates 

The Airy function Ai(x) and its derivative have along the positive real 
axis, and more generally as |2;| — )• oo with | arg z\ < tt — 5 for any fixed 5 > 0, 
asymptotic expansions given in full in [1, 10.4.59 and 10.4.61]; we need only 
the leading terms: 

Ai(z) ~ ^^z~y'e-h"\ (A.l) 
Ai'(z)~-^zVVi.3/^. (A.2) 

In particular, along the imaginary axis, for y € M, 

\K\{^y)\r.^\y\-y^e'i\y^''\ (A.3) 

Along the negative real axis, Ai and Ai' oscillate and have zeros; by fH 
(10.4.60), (10.4.62)] (where further terms are given), we have the asymptotic 
formulas 

Ai(-x) = ^-V2a;-i/4(sin(2^3/2 ^ ^ ^ 0{x'^/^), (A.4) 
Ai'(-x) = -7r-i/2xi/4(cos(fa;3/2 + |) + o(l)). (A.5) 
and more generally, as |2;| — t- oo in any domain | argz| < |7r — (5, 

M{-Z) = ^-l/2z-l/4(^sin(^z3/2 + ^) (1 + 0(|z|-3/2)) ^ 0{\z\-^l^)^ , 

(A.6) 
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2 



+ 9(l + 0(|zr3/2))+o(|.|-=^/2)). 

(A.7) 



For the companion Airy function Bi(x), we use the following estimate [U 
10.4.63], valid along the positive real axis, and more generally as \z\ — >• oo 
with I argzj < 7r/3 — 5, 



BiU)~7r-V2^-i/4ei^'''. 



(A. 



(A.9) 



On the negative axis, we have [1, 10.4.64], 

Bi(-x) = 7r-i/2x-i/4(cos(|x3/2 + |)+o(l)). 

Lemma A.l. For every 5 > there exists c = c{5) such that if \ argz| < 
-n — 5 and x > 0, then 

Ai{z) 



Proof. We may assume that < 5 < tt/2. It follows from (A.l) that for 
I arg z\ < TT — 6, 

\Ai{z)\ X (l + |z|)-i/^exp(-|Re(z3/2)). (A.IO) 

Since | arg(2; + x)\ < \ arg(2;)| < n — S, ( A.lOp can be used for z + x too. If 
Rez < 0, then |z+x| > | Imz| > (sin 5)1^1, and if Rez > 0, then |z+x| > \z\. 
Hence, ( A.10| ) implies 

Ai(z + x) 



Ai(z) 



|Rc{(2+x)3/2)+|Re(^3/2) 



Further, 



Re{{z + x f/^) - \ Re(z=^/2) ^ ^ ^)i/2 

3 Jo 



> cos ■ 



t|V2dt. 



Moreover, |-2 + t| > ci(|z| + t) by elementary geometry (use, e.g., the sine 
theorem on the triangle with vertices 0, z, —t); hence j^+tj-'^/^ > C2(|2:|+i)^/^, 
and the result follows. □ 

Recah that Rn := (|7riV)^/l 

Lemma A. 2. Let x > be fixed. Assume that z = Rj^e^^ with 9 £ [f , ^] 
and N > 1. Then (with implicit constants depending on x but not on N or 
0), 



and 



\M{z) 



Ai{z + x) ^^(^-.Rl/^\e-^\/^ 
Ai{z) 



|^|-l/4g|lmf(-^)3/2| ^ ^-l/4g|ij3/2sm{3|^|/2)^ 



(A.ll) 
(A.12) 
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Proof. By symmetry, it suffices to consider 9 = n + with < 99 < 7r/2. 
Thus z = —Rnc^'^. We have, by Taylor's formula. 



and thus 



Im - {—z — X 

—3 /2 

If 99 > Rj^ , this shows that, for large N ^ 



(A.13) 
(A.14) 



> 0.2, 



and thus, by (A.6), 



|Ai(z + x)\ X |z + x| 



-1/4 



sm 



-z - xf/^ + 



xii^^/%xp(lm|(-z-x)3/2). (A.15) 

In particular, taking x = we obtain (A. 12) for these ip. Moreover, com- 
paring ( |A.15 ) with the special case x = 0, and using (A.14), 

Ai(2: + x) 



Ai(z) 



exp(lm|(-z - x)^/^ -Im|(- 



^3/2 



( dV2 ■ 'P\ 



This yields ( |A.llD for i?^^''^ < 9? < vr/2, since sin(99/2) > (^/vr = |6i - vrl/vr. 

For < < R"^'"^, we find from ^KI^ Im|(-z-x)3/2 = 0(1), so 
yields Ai(z + x) = 0{\z\~^/'^). Similarly, by our choice of Rn, 
'2 



Re(|(-.)^/^) = l^/'cosf = '^rT + 0{R,"')=.N + 0{R,'"). 

Hence, for large at least, |sin(|(— z)'^/^ + |)| x 1, and, by ( |A.6[ ) again, 

Ai(z) = e){\z\-^/'^)\ consequently Ai(2; + x)/Ai(z) = 0(1). These are the 

1/2 3/2 
results we want in this case because Rn <~p = 0(1) and i?„ 99 = 0(1). □ 

Recah also the entire functions Gi and Hi defined by P (10.4.44), (10.4.46)]: 
Hi(z) := vr"^ y exp(-^t^ + ztj dt 
Gi{z) := Bi(z) - Hi(z). 



(A.16) 
(A.17) 

The integral defining Hi evidently converges for all complex z. There is a 
similar integral yielding Gi fl\ (10.4.42)]: 

1 



Gi(z) 



TT 



sin( ^t^ + zt] dt, 



z G 



(A.18) 



but this integral converges only conditionally for real z, and not at all for 
z^R. 

We have the asymptotics, as x — )• 00, [U (10.4.91)] 

Hi(-x) ~ vr-^x-^ (A.19) 



26 SVANTE JANSON, GUY LOUCHARD, AND ANDERS MARTIN-LOF 

More generally, for Re 2: < 0, expanding exp(— in (A. 16) yields, for any 
L > 0, 



Hi(z) 



-vr 



L-l 



3^ 



+ 0(|Rez|^3L-i^ _ 



In particular, for any 5 > and argz G (f + 5, 



f -5), 



Hi(z) = 0(|Rez|-^^ 
More precisely, for such z, 



0(|; 



1-1^ 



Hi(z) 



-vr 



-iz-i + o(|z|-4) 



We further define, using Ai{x)dx = 1/3 [1, 10.4.82], 



AI(z) 



+00 2 



Ai(t) dt 



Ai(t) dt; 



(A.20) 

(A.21) 
(A.22) 

(A.23) 



this is well-defined for all complex z and yields an entire function provided 
the first integral is taken along, for example, a path that eventually follows 
the positive real axis to +00. Note that Al'(2:) = —Ai{z) and that AI(0) = 
1/3. Along the real axis we have the limits, by [1, 10.4.82-83], 

lim AI(x) = 0, (A.24) 



/oo 
Ai(x) dx = 1. 
-00 



(A.25) 



In terms of the functions Gi and Hi, we have, see [H 10.4.47-48], 

Al{z) = 7r(Ai(z)Gi'(z) - Ai'{z)Gi{z)) (A.26) 

= 1 + 7r(Ai'(z)Hi(z) - Ai(z)Hi'(z)) . (A.27) 

We have, see [HI Appendix A] and (for real z) [I, 10.4.82-83], as \z\ — 00 
with I arg(2;)| < it — 5, 



Al(z) 



2^AF 



(A.28) 



and, along the negative real axis, and more generally for —z with | arg(2;)| < 
27r/3-<5, 

AI(-z) = l-7r-V2^-3/4j^^^gJ^2^3/2^^^(i^o(|^|-3/2))^0(|^|-3/2)^_ 

(A.29) 

We have the following estimates. 
Lemma A. 3. (i) For every fixed 6 > 0, if \ argz| < ir — 6 and \z\ > 1, then 



Ai^ 
Ai{z) 
AI(z) 
Ai{z) 



,1/2 



-1/2 
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;ii) For z = Rnc''^ with 9 G [f , f] and N >1, 
Ai'(z 



Ai(z) 
Aliz] 



Ai{z) 



O (|z|-i/2 + |Ai(z)|-i) = O + |z|V4e-§|lm{-.)3/2|^| 



Proof. Part (i) follows from (A.l), (A.2) and (A. 28). 



For (ii), we use (A. 12), (A. 7) and (A.29). 



□ 



Let afc, k > 1, denote the zeros of the Airy function; these are all real and 
negative, so we have > ai > 02 > By [U (10.4.94), (10.4.105)], 



^-'^^^-'^ ^'(1 + 0(k-')) ~ )"' . k'lK (A.30) 



Thus, by (|A^, see also P (10.4.96)], \A^ , ( |Al9l ) and \A?n\ , 

|Ai'(a,)|x|afc|i/^xfci/6^ (A.31) 

|Bi(a,.)| ^ lafcl-i/^ X k-^'\ (A.32) 

|Hi(a,.)| lafcl-i X A;-2/3. (A.33) 

|Gi(afc)| X |afc|-i/4 X A;-^/^, (A.34) 



Appendix B. Some Airy integrals 

Integrals of the Airy functions (and their derivatives) times powers of x 
are easily reduced using the relations Ai"(x) = xAi(3;) and Bi"(x) = xBi(x) 
and integration by parts. We have, for example, using also the definition 



(A.23) 



Ai(x) dx 
xAi(x) dx 



AI(x) 

Ai"(x) dx = Ai'(x) 



(B.l) 
(B.2) 



/ x^Ai(x) dx = / xAi"(x) dx = xA.'(x) - / A.'(x) dx = xAi'(x) - Ai(x) 

(B.3) 

and in general the recursion 

j x"Ai(x) dx = y x"-iAi"(x) dx = x"-iAi'(x) - (n - 1) j x"-2Ai'(x) dx 
= x"-^Ai'(x) - (n - l)x"-2Ai(x) + (n - l)(n -2) j x"-3Ai(x) dx. (B.4) 
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Integrals of products of two Airy functions and powers of x can be treated 
similarly, see ; we quote the following (that are easily verified by differen- 
tiation): 

j Ai{xf dx = xAi{xf - {Ai'{x)f, (B.5) 

j xAiixf dx = l {x'^Ai{xf - x{Ai'{x)f + Ai'(x)Ai(x)) (B.6) 

J x^Ai{xfdx = l{x^A[{xf - x'^{Ai{x)f + 2xAi'(x)Ai(x) - A[{x)^) 

(B.7) 

and in general the recursion 

j x"Ai(x)2 dx = ^ (^x"+iAi(x)2 - x"(Ai'(x))2 + nx''-^A[{x)Ai'{x) 

n(n — 1) „_2,../ n(n — l)(n — 2) f „_o,,., ,,2 , A /t^ 
^— — ^x" ^(Ai(x))^ + ^ ^ / ^ ^(Ai(x))^dx j. (B.8) 

By the same method, we can also treat products involving two different 
translates of Airy functions; this gives for example the following (again, 
these are easily verified by differentiation), if a 7^ 6 and c = (a + b)/2: 

j Ai{x + a)Ai(x + b)dx = {Ai{x + a)Ai(x + b) - Ai{x + a)Ai'(2; + 6)) , 

(B.9) 



y"(x + c)Ai(x + a)Ai(x + b) dx 



^ ' 'a-b){x + c){Ai'{x + a)Ai{x + b) - Ai(x + a)Ai'(x + 6)) 



(a - 6)2 

- 2{x + c)Ai(x + a)Ai(x + b) + 2Ai{x + a)Ai'(x + b) 
2 

+ 7 (Ai'(x + a)Ai(x + 6) - Ai(x + a)Ai'(x + b)) , 

(a — b)-^ 

(B.IO) 

and the recursion 



(x + c)"Ai(x + a)Ai(x + 6) dx = 
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/ ,N9 1 (« - b){x + c)"'(Ai(x + a)AUx + b) - Ai(x + a)Ai(x + b)) 
(a — o)^ \ 

- 2n{x + c)"Ai(x + a)Ai(x + 6) + 2n{x + c)""^Ai'(a; + a)Ai'(x + b) 

- n{n - l)(x + c)"^^(Ai'(x + a)Ai(a; + b)+ Ai{x + a)Ai'(2; + b)) 
+ n(n - l)(n - 2)(x + c)""^Ai(x + a)Ai(x + b) 

+ 2n(2n - 1) j {x + c)""^Ai(x + a)Ai(x + b) dx 

- n{n - l)(n - 2)(n - 3) j {x + c)""'^Ai(x + a)Ai(x + b) dx^ . 

(B.ll) 

In particular, if ak and are zeros of Ai, with k ^ the formulas above 
yield, recalling the rapid decay ( A.l[ ) and (A. 2) at oo: 



rOO /'OO 

/ Ai(x + Ofc) dx = Ai(x) dx = Al(ajfc) = — 7rAi'(afc)Gi(afc), 

Jo Jak 

l-OO 

/ xAi(x) dx = — Ai'(afe), 

/>oo 

/ x^Ai(x) dx = — afcAi'(afc), 
j xAi(x + afc) dx = -Ai'(afc) - afcAI(ofc), 

JO 



Ai(x + Ofc)^ dx 



oo 



oo 



Ai(x)2dx = (Ai'(ofc))2, 
xAi(x)^dx = iafc(Ai'(afc))^, 



^Ai(x)2dx = i4(Ai'(afc))2, 



;>oo 

/ xAi(x + Ofc)^ dx = — 

JO 

/•oo 

/ Ai(x + afc)Ai(x + a^) dx = 0, 

JO 
/•oo 

/ xAi(x + afc)Ai(x + a^) dx 

JO 



|afc(Ai'(afc))' 



{ak - ae)' 



-,Ai{ak)A\'{a(). 



(B.12) 
(B.13) 
(B.14) 
(B.15) 
(B.16) 
(B.17) 
(B.18) 
(B.19) 
(B.20) 
(B.21) 



More generally, for arbitrary complex a 7^ 6, by ( B.9| ) and (B.IO), again 
using ( 1X21 ), 

/•oo 1 

/ Ai(x + a)Ai(x + 6) dx = -(Ai(a)Ai'(6) - Ai'(a)Ai(6)), (B.22) 

Jo a-b^ 

f°° a + b 
/ xAi(x + a)Ai(x + 6) dx = -^Ai(a)Ai(5) 

Jo (« - by 



(a -by 



,Ai'(a)Ai'(6) 
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+ 



(a - by 



■(Ai(a)Ai'(6) - Ai'(a)Ai(6)). 



In particular, 



Ai(x)Ai(x ■ 
xAi(a;)Ai(x ■ 



ak) dx 
cLk) dx 



Ai(0)Ai'(afc) 



-ak 



(B.23) 



(B.24) 



-^Ai'(0)Ai'(afc) - ^Ai(0)Ai'(afc). (B.25) 

Remark B.l. The Sturm-Liouville operator Tf{x) = —f"{x) + xf{x) on 
[0, oo) with the boundary condition /(O) = has the eigenfunctions Ai(x + 
ak) with eigenvalues —ak = \ak\- (B.20) thus expresses the orthogonality 
of the eigenfunctions which also follows by general operator theory; in fact, 
the operator T is self-adjoint and these eigenfunctions form an orthogonal 
basis in L^(0,oo). The corresponding ON basis is by (B.16) given by the 
functions Ai(a; + afc)/Ai'(afc), k > 1. 

For example, (2.10) is the expansion of G{2~^^^x) in this ON basis, with 
coefficients 7rHi(ajt), which yields (4.5) by Parseval's formula. Similarly, 
Ai{x) has the expansion (convergent in L^) 



Ai(x) = J2 



Ai(0) 



-Ai(x + ak), 



^ \ak\Ai'{ak) 

where the coefficients are given by ( |B.24 ), and thus, using (B.5), 



AiW 
Ai(0)2 



/Ai(x)Y 
VAi(O); 



dx 



oo ^ 

k=l 



(B.26) 



(B.27) 



We will also use the Laplace transform of the Airy function, which is 
easily found. (Taking z imaginary, we obtain the Fourier transform e*^ 
of Ai; this is sometimes taken as the definition of Ai, see e.g. |10^ Definition 
7.6.8].) 



Lemma B.2. If Kez > 0, then 

/oo 
e^*Ai(t) dt 
-oo 



Proof. By (A.l) and (A.4), the integral converges absolutely for every z 
with Kez > 0, and thus the integral is an analytic function of z in the right 
halfplane, say F{z). We have, for Rez > 0, by Ai"(t) = tAi(t) and two 
integrations by parts. 



F'{z) 



/OO /"OO 
te^*Ai{t)dt= / e^*Ai"(t)dt 
-OO J — OO 

/OO rOQ 
ze^*Ai'(t) dt = / z'^e^*A\{t)dt = z'^F{z). 
-oo J —oo 



Hence, F{z) = Ce^^^^ for some C. 
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/oo POO POO POO POO 

e^*Ai(t)dt= / ze'' A[{t)dtds= e" / Ai(t)dtdu. 
-oo J —oo Js J —oo Ju/z 

(B.28) 



oo 



Since AI(x) := r^Ai(i) dt — t- as x — t- +oo and AI(x) — t- 1 as x - 
by by ( A.24| )-( A.25 ), dominated convergence shows that, letting z \ in 

/oo pO 
e"l[u <0]du= e^'du = 1. 
-oo J —oo 



□ 



Remark B.3. By Fourier inversion we find, for any a > 0, 



Ai(t) 



1 

27ri 



(B.29) 



By analytic extension, this holds for any complex t. 

Appendix C. An integral equation for fr{t) 

We give here another approach, based on Daniels [personal communica- 
tion, 1993], to find the density function of the defect stopping time r = t^, 
which was the basis of our development in Section [3j Unfortunately, we have 
not succeeded to make this approach rigorous, but we find it intriguing that 
it nevertheless yields the right result, so we present it here as an inspiration 
for further research. 

Let as in (3.1) be the first passage time of W{t) to the barrier b{t) := 
—X — t^/2, where x > is fixed. Let frit) be the (defect) density of r, and 



/\/27rt the density of W{t). 



The first entrance decomposition of W{t) to the region w < b{t) gives the 
integral equation (using the strong Markov property) 

(j){w;t) = / fT-{u)4>{w — b{u);t — u) du 
Jo 

for w < b{t). Letting w b{t) we get the equation 



mt)-t)= / ^(n)0(6(O-6(n);t-n)dn 



(C.l) 



for t > 0. (Similar arguments using the last exit decomposition, which leads 
to another functional equation involving also another unknown function, are 
used by Daniels [6] and Daniels and Skyrme [8].) Since 

b{t) - b{u) = (n^ - t2)/2 = -(t - u){t + u)/2, 



we have by (C.l ) 



^-{t-u){t+u)^/8 
friu) . - du. 



V^27r(t - u) 



(C.2) 
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The exponents can be written as —{t — u){t + m)^/8 = + + {t — 

uf/2A and -{x + t^/2f/2t = -x'^ /2t - xt/2 + t^/2A - 4^6, so the integral 
equation ( C.2[ ) can be transformed into 



x'^ /2t-xt/2+t^ /2i ft Jt-uy^/2A 

' fr{u)e'' /^-^ du. (C.3) 



V2M Jo \l2Ti{t - u) 

This is a convolution equation of the form 

h{t) = / g{u)k{t — u)du 
Jo 

with 

g{t) := ^(^)e*'/^ 

x'^/2t-xt/2+t^/24: -b{t)^ /2t+t^ /6 



e ' ' ' ' e 

h{t) : = 



V2TTt V2TTt 

pt^/24: 

kit) , 



s/2-Kt 

If the Laplace transforms ^(s) := e~^*g{t) dt etc. were finite for Res 
large enough, we could get the solution from g{s) = h{s)/k{s). However, the 
factors e* /^"^ in h{t) and k{t) grow too fast, so h{s) and k{s) are not finite for 
any s > and this method does not work. Nevertheless, if we instead define 
his) and kis) by his) := /;+~' e-*/i(i) di and fc(s) := /;+^; e-^*fc(t) di, 
integrating along vertical lines in the complex plane with real part a > 0, 
then the formula gis) = his) /kis) yields the correct formula for gis) and 
thus for git). (Note that h and k can be seen as Fourier transform of h and 
k restricted to vertical lines. The value of cr > is arbitrary and does not 
affect h and k.) Let us show this remarkable fact by calculating /i(s) and 
kis). 

Consider first and express the Gaussian factor by Fourier inversion: 

V2TTt Ja-ooi 27rZ 



The exponent 6(t)u + tu^/2 can then be written as —ux + u^/Q + it — u)'^/6- 
t^/6, so that 

(J+OOj SI 

ux+u^/6+(t-u)'-y6 



a—ooi 
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and, choosing ai > a > and letting a2 ■ 

f(Tl+OOj 



0"! 



cr, 



h{s) 



a+ooi 



— ooi 
cr+ooi 



st—ux+u^/6+{t—u)^/6 



ooi 

(T2 + OOi 



-s{u+v)—ux+u^ /6 



du dt 
~27d~ 

dv du 



0"2— 002 
(T2 + OOi 



(72—001 



^-sv+v^ /6^-{s+x)u+u^ /6 



2-Ki 

dv du 
2m 



= 2mK\[c{s + x))Ai(cs)c^, 

with c := 2^3^ using ( |B.29| ). 

Since k{t) is obtained by putting x = in we get directly 

k{s) = 2'KiA\{csfc^ 

and hence 

h{s)/k{s) = ki{2^f'\s + x))/Ai{2^''h). 

This is indeed the Laplace transform of g{t) = /r(i)e*^^^ given in (3.8), 
which by inversion yields the formulas (3.7) and (3.6) for frit)- 

It seems likely that it should be possible to verify the crucial formula 
g{s)k{s) = h{s) by suitable manipulations of integrals, which would give 
another proof of the formulas (3.6)-(3.8) for frit) and g{t). For example, if 
we define, for Ret > 0, 

/■Ret 

F{t):=h{t)- g{u)k{t-u)du 







(note that F is not analytic), then F{t) = for real t > 0, and it is easily veri- 
fied that the equation g(s)k{s) = h(s) is equivalent to j^^^l e~^^F(t) dt 
as cj — )• oo. However, we do not know how to verify this directly. We there- 
fore leave finding a direct proof of g{s)k{s) = h{s) as an open problem. 



Appendix D. An alternative derivation of (3.8) 



A proof of the formula (3.8) for the Laplace transform of the density of the 



passage time r is given by Groeneboom f9j (in a more general form, allowing 
an arbitrary starting point and not just t = 0, or, equivalently, linear term bt 



in (1.1) or (3.1); for simplicity we do not consider this extension). His proof 
uses partly quite technical methods. For the service of the reader we here 
present an alternative proof based on the same ideas but from a different 
point of view; we believe that this yields a more straightforward proof for 



our purposes. As discussed in Remark 3.3, this implies (3.6) and (3.4), so it 



3.2 



gives us a self-contained proof of the central Lemma ^ 

As in [ini [IS] we consider the process with drift t^/2 defined by X{t) : 
X + 1^/2 + W{t), so that T defined by (with /3 = 1/2 as in (|0|-([3^) 
is the first hitting time ol X{t) = 0. If Px is the probability measure (on 
the space C[0,oo)) corresponding to X{-), and Qx that corresponding to 
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X + W{-), then the Cameron-Martin formula tells us that, considering the 
restriction to a finite time interval [0, t] , the Radon-Nikodym derivative is, 
using dX{s) = sds + dW{s), 



dQx{X) """"PV 2 JO 



1 /•* f{dX{s) - sdsf dX{s) 



ds 

1 



ds 



exp I / sdX(s) s ds 

\Jo 2 Jq 

exp (-^ + tX{t) - X{s) ds^ 



using integration by parts. Hence, letting denote expectation with re- 
spect to the Wiener measure Qx, 



P.(r>t)=E.(^;r>t) =e-*='/6E.(e*Wo'^Wd-^>t) (D.l) 



dQx 



(cf. jl6f Theorem 2.1] and [9, Lemma 2.1]). We introduce the "Green func- 
tion" 

F(x,y,t) :=E, (e-/o - y); r > t) , 

where 5 is the Dirac delta function (formally, F{x,-,t) is defined as the 



density of the corresponding occupation measure); (D.l) then says that 



P,(t >t) = e-*'/6 / F{x, y, t)e'y dy. (D.2) 
Jo 



The Feynman-Kac formula tells us that for fixed y > 0, F{x, y, t) is the 
fundamental solution of the equation 



dF _ 1 d'^F 
~dt ~ 2 5x2" 



xF, x,t > 0, 



(D.3) 



with the boundary conditions F{0,y,t) = F{oo,y,t) = 0. By time-reversal 



(or by symmetry of the resolvent R{x,y,z) in (D.7) below), F{x,y,t) 
F{y,x,t)] hence also 



dF 1 a^F 



(D.4) 
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Differentiating (D.2) under the integral sign and using the Feynman-Kac 
equation (D.4), we obtain (cf. O Lemma 2.2]) 



F{x,y,t)e'ydy 



oo 





1 
2 

1 dF 
2d^ 



^F(x, y, t) - _ t) ) e'y dy 



te^yFix,y,t)-e'y^^^^^ 



{x,0,t) 



(D.5) 



since F{x,0,t) = and F and its derivatives decrease rapidly as y — )• oo. 

The Laplace transform of F{x,y,-) is the resolvent, defined at least for 
Rez > 0, 



POO 

Rix,y,z) := / e-^*F(x, y, t) dt. 
Jo 



(D.6) 



The Feynman-Kac equation is a second order differential equation, and the 
theory of such equations (see also [HI Appendix C; in particular (379)] for 
this particular case) tells us that 



Rix,y,z) 



^ipo{x;z)ipoo{y;z), 0<x<y, 
^ipo{y;z)ipoo{x;z), < y < x, 



(D.7) 



where ^'0(2;) = ^o{x; z) and ipoo{x) = ipoo{x; z) are solutions of the differen- 
tial equation 



(x) — Xip{x) = z(p{x) 



with boundary conditions V'o(O) = 0, 9300(00) = 0, and w is the Wronskian 
w := ifoo{x)ip'Q{x) — LpQ{x)ip'^{x) (which is constant in x). Note that since 
(^o(O) = 0, 



W = (/?oo(0)(/?o(0). 



(D.9) 



The differential equation (D.8) has two linearly independent solutions 
A{x + z) and B{x + z) with A{x) := Ai(cx) and B{x) := Bi(cx) with 
c := 2^/^, and in terms of these we have (up to arbitrary constant factors) 



Lpo{x] z) = A{z)B{x + z)- B{z)A{x + z), 
LPoo{x] z) = A{x + z). 



(D.IO) 
(D.ll) 
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We integrate (D.5) (multiplied by e ) and differentiate (D.6) and obtain, 



using also ( |D.7D , ( |D.9D and ( |D.llD 



^a;,0,t)dt = - — (x,0,2) 



which is (3.8 ). 



dy 



2dy 



A{x + z) 
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